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ABSTRACT 

We show that the coronal heating and the fast solar wind acceleration in the coronal holes are natural 
consequence of the footpoint fluctuations of the magnetic fields at the photosphere, by performing one- 
dimensional magnetohydrodynamical simulation with radiative cooling and thermal conduction. We 
initially set up a static open flux tube with temperature 10 4 K rooted at the photosphere. We impose 
transverse photospheric motions corresponding to the granulations with velocity (dv±) = 0.7km/s and 
period between 20 seconds and 30 minutes, which generate outgoing Alfven waves. We self-consistently 
treat these waves and the plasma heating. After attenuation in the chromosphere by ~ 85% of the 
initial energy flux, the outgoing Alfven waves enter the corona and contribute to the heating and 
acceleration of the plasma mainly by the nonlinear generation of the compressive waves and shocks. 
Our result clearly shows that the initial cool and static atmosphere is naturally heated up to 10 6 K 
and accelerated to ~ 800km/s. 
The mpeg movie for fig.l is available at 

\protect\vrule widthOpt\protect\href {http : //www-tap . scphys .kyoto-u. ac . jp\string~stakeru/research/ suzuki_200 

Subject headings: magnetic fields - plasma - magnetohydrodynamics - Sun : corona - solar wind - 
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1. INTRODUCTION 

A certain portion of the solar surface is covered by the 
coronal holes. Their characteristic features are that the 
open magnetic fields are emanating into the interplan- 
etary space and the hot plasma is streaming out. The 
high speed solar winds (800km/s at ~ 1AU) are accel- 
erated from the polar coronal holes whic h stably exists 
excep t during the solar maximum phase ((Phillips et al.l 
Il995ft The source of the energy for heating and accelerat- 
ing the plasma is believed to be in the surface convection; 
this energy is lifted up through the magnetic fields and its 
dissipation results in the plasma heating. In this sense, 
the problem of the coronal heating and the solar wind 
acceleration is to solve how the solar atmosphere reacts 
to the footpoint motions of the magnetic fields. 

In the coronal holes the Alfven waves can play an im- 
portant role because it can travel a long distance to 
both heat the corona and accelerate the solar wind. 
The Alfven waves are excited by steady transverse mo- 
tions of the field lines at the photosphere (e.g. Cran- 
mer & van Ballegooijen 2005), while they can also be 
produced by reconnection e vents above the photosphere 
((Axford fe Mckenzie! 11997(1 . Although the latter pro- 
cess might be responsible for the generation of the high- 
frequency (up to 10 4 Hz) ioncyclo tron wave h ighlighted 
for the heating of the heavy ions llKo hl et al. 1ll998f) . it 
is difficult to heat the protons <(CraninerTl2000() ~]3e- 
cause we focus on the heating of the major component 
of the plasma, we study the effect of the low-frequency 
(< 0.1Hz) Alfven waves generated by the steady foot- 
point fluctuations. 

Although the Alfven waves have been intensively in- 
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vestigated, previous studies have a couple of limitations. 
Most of the calculations consider the waves from the fixed 
'coronal base' of which density is ~ 9 orders of magni- 
tudes lower than that of the photosphere (e.g. Oughton 
et al.2001; Ofman 2004). However, the corona dynam- 
ically interact with the chrom osphere by chromospheric 
evaporation iH ammer | (19821) and wave transmission 
l(Kudoh fc Shibata Ill999|) so that the coronal base varies 
time-dependently. Some calculations including the chro- 
mosphere even adopt an ad hoc heating function (e.g. 
Lie-Svendsen 2002). In contrast, we self-consistently 
treat the transfer of the mass/momentum/energy by 
solving the wave propagation from the photosphere to 
the interplanetary region. Our aim is to answer the prob- 
lem of the heating and a cceleration in the coronal hole s 
by the forward approach ((Cudiksen fc Nordlund 1 12005) . 

2. METHOD 

We consider one-dimensional (1-D) super-radially open 
flux tube, measured by heliocentric distance, r. The com- 
putation domain is from the photosphere (r = li?s) with 
density, p — 10~ 7 g cm -3 , to 65i?s (0.3AU), where i?s is 
solar radius. Radial field strength, B ri is given by con- 
servation of magnetic flux as 



B r r f(r) = const. 



(1) 
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where we adopt the same function as in lKopp fc Holzerl 
l(1976() for superradial expansion, f(r), but consider two- 
step expansions: the flux tube initially expands by a fac- 
tor o f 30 at r ~ l.O lifc corresponding to 'funnel' struc- 
ture l(Tu et all2005|) . and followed by 2.5 times expansion 
at r ~ 1.2i?s due to the large scale magnetic fields. B r at 
the photosphere is set to be 161G. These give B r = 5G 
at r = 1.02.Rs (low corona), and B r = 2.lG(R s /r) 2 in 
r > 1.5i?s- We give transverse fluctuations of the field 
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Fig. 1. — Time evolution of the atmosphere. From top to bottom, 
outflow speed, t) r (km/s), temperature, T(K), density, p (g cm -3 ), 
and rms amplitude of transverse velocity, (dv±)(km/s) are plotted. 
Dashed, green, blue, and red solid lines are results at t = 0, 20, 340, 
& 2573 mins., respectively, whereas results at t = do not appear 
in v r and (dvj_) because both equal to 0. The results are averaged 
by integrating for 3 minutes to take into account the effects of the 
exposure time for comparison with observations. 



line by the granulations at the photosphere, which ex- 
cite Alfven waves. We consider the fluctuations with 
power spectrum, P(y) oc v^ 1 , in frequency between 
6 x f CT 4 < v < 0.05Hz (period of 20 seconds — 30 min- 
utes), and root mean squared (rms) average amplitude 
(dvj_) — 0.7km/s corresponding to observed velocity am - 
plitude ~ lkm/s l|Holweger. Gehlsen. fc Ruland II197S|) . 
At the outer boundary, outgoing condition is imposed for 
all the MHD waves (Suzuki & Inutsuka 2005) to avoid 
the unphysical wave reflection. 

We treat the dynamical evolution of the waves and the 
plasma by solving ideal MHD equations with the relevant 
physical processes : 



d P , P d 1 2 1 \ n 



dt 



2 f dr 



(2) 



dv r 
} ~dt 



dp 
dr 



d 



8-7rr 2 / dr 



(r 2 !Bl) 



2r 2 f dr 



B r d 



GM S 



-( — — 



47r dr 
GMq 1 d 
r r 2 f dr 



y/fB±) 



(3) 
(4) 



B 2 



8tt 



8B± 

~dT 



i a 

r^/J dr 



[ry/f(v±B r - v r B±)], 



(5) 



(6) 



f 



where p, v, p, e, B are density, velocity, pressure, specific 
energy, and magnetic field strength, respectively, and 
subscript r and _L denote radial and tangential compo- 
nents. ^ and denote Lagrangian and Eulerian deriva- 
tives, respectively. G and Ms are the gravitational con- 
stant and the solar mass. F c is thermal conductive flux 
and <7r is radiative cooling llLandini &: Monsignori-Foss 



112221 lAnderson fc Athavll1989t iMorivaan et a,Ul200 
Note that the curvature effects appear as ry/J terms 
instead of r for the usual spherical coordinate. We adopt 
2nd-order MHD-Godunov-MOCCT scheme for updating 
the physical quantities (Sano & Inutsuka 2005), of which 
an advantage is that no artificial viscosity is required 
even for strong MHD shocks. We fix 14000 grid points 
with variable sizes in a way to resolve Alfven waves with 
v = 0.05Hz by at least 10 grids per wavelength to reduce 
the numerical damping. 

The advantages of our calculations are automatic 
treatments of the wave propagation/dissipation and the 
plasma heating/acceleration without an ad hoc heat- 
ing function but with the minimal assumption over the 
broadest regions with the huge density contrast. We em- 
phasize that this is one of the most self-consistent sim- 
ulations for the solar wind acceleration at present. The 
disadvantage is the 1-D MHD approximation, which will 
be discussed later. 

3. RESULTS AND DISCUSSIONS 

We start the calculation from the static and cool at- 
mosphere with temperature, T = 10 4 K; we do not set up 
the corona and the solar wind at the beginning. Figure 
1 (on-line movie available) shows how the coronal heat- 
ing and the solar wind acceleration are realized by the 
low-frequency Alfven waves. We plot u r (km/s), T(K), 
p(g/cm 3 ), and (dv±)(km/s) averaged from v± as a func- 
tion of (r - R s )/Rs at different time, t = 0,20,340 
& 2573 minutes. As time goes on, the atmosphere is 
heated and accelerated effectively by the dissipation of 
the Alfven waves. Temperature rises rapidly in the inner 
region even at t = 20 minutes, and the outer region is 
eventually heated up by both outward thermal conduc- 
tion and wave dissipation. Once the plasma is heated 
up to the coronal temperature, mass is supplied to the 
corona mainly by the chromospheric evaporation due to 
the downward thermal conduction. This is seen in tem- 
perature structure as an inward shift of the transition re- 
gion, which is finally located around r — Rg = 6 x 10 _3 i?s 
(~4000 km). As a result, the coronal density increases 
by two orders of magnitude. While the wind velocity 
exceeds 1 000km/ s at t = 340 minutes on account of the 
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Fig. 2. — Comparison of the result at t = 2573 mins (red solid 
lines) with observations summarized below. Quantities in four 
panels are the same as in Figure 1 except third panel showing 
electron density, N e (cm~ 3 ), instead of p(g cm -3 ). Scales in both 
horizontal and vertical axi ses are changed fro m fig.l. a: Vertical 
error bars with green x's ITeriaca et al. 2003) and blue triangles 
IZangri lli et al. 2002) arc proton outflow speed in polar regions by 
SoHO. Blue square with errors is velocity by IPS measurements 
avera ged in 0.13 - 0.3AU of high-latitude regionsl Koiima ct al. 
120041) . Crossed bars with and witho ut circles a re measurements 
by VLBI and IPS(EISCAT) IGrall et al J 1199(1 ■ Vertical error 
bars wi th circles are data based on observation by SPARTAN 
201-01 IHabbal et aH[l995l). b: Pink circles arc electron temper- 
ature by CDS /SoHO IFludra. Del Zanna. fc BromagH I199SD . c: 
Pink circles and [blu e] stars are observatio ns by SU MER/SoHO 
IWilhelm et alj 119911) and by CDS/SoHO JTeriaca et alj|2003T) . 
respectively. [Gre enl triangles QTeriaca"1itHlrtTl27)o'3T) "lajid squares 
ILamv et al.l 119971) are observations by LASCO/SoHO. d: Blue 
circles are non -therma l broadening inferred from SUMER/SoHO 
measurements (Bancricc ct al. 1998). Cross hatched region is em- 
pirically der ived non-thermal broadening based on UVCS/SoHO 
observation (Esser et al. 1993). Green error bars with circles are 
transvers e velocity fluctuations derived from IPS measurements by 
EISCAT (Canals et al. 200^). 



initial low density, it gradually settles down to < 1000 
km/s as the density increases. (dv±) also settles down to 
the reasonable value at the final stage. 

We have found that the plasma is steadily heated up 
to 10 6 K in the corona and flows out as a transonic wind 
with v r ~ 800km/s at the outer boundary (=0.3AU) 
when the quasi steady-state behaviors are achieved after 
t > 1800 minutes. This is the first numerical simula- 
tion which directly shows that the heated plasma ac- 
tually flows out as the transonic wind, from the ini- 
tially static and cool atmosphere, by the effects of the 
Alfven waves. The sonic point where v r exceeds the lo- 
cal sound speed, c s = y^5p/3p, is located at r ~ 2.5i?s 
and the Alfven point for the Alfven speed, B r / 
is at r ~ 24i?s- The obtained proton flux at 0.3AU 
(2 ± 0.5) x 10 9 (cm _2 s~ 1 ), corresponding to 
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100 



(r-Rs)/Rs 



Fig. 3. — Top : S c (see text) of outgoing Alfven mode (red), in- 
coming Alfven mode (blue), and outgoing MHD slow mode (green) 
at t = 2573 mins. Hatched region indicates the chromosphere and 
low transition region with T < 4xl0 4 K. Bottom: Normalized wave 
amplitude, (<5«a+)/ c S (solid), (Sva +)/va (dashed [red solid]), and 
((Jug +)/cs (dotted [green solid]). 



In fig. 2 we compare the result at t = 2573 minutes 
with recent observation in the high-speed solar winds 
from the polar regions by Solar fe Heliospheric Obser- 
vatory (SoHO) iZangrilli et all2002t ITeriaca et alJ |2003: 



Fludra. Del Zanna. fe Bromagdll99S : ILamv et al.Ml997 



Wilhelm et all 119981: iBaneriee et all (l998t lEsser et alJ 
1999) and Interplanet a ry Sc in tillation (IPS) mea- 
surements llOrall et all (l99ft iKoiima et al~l 12004 



iCanals et alJl2002j) . The fi gure shows that our forward- 
approach simulation naturally form the corona and the 
high-speed solar wind which are observed by the dissi- 
pation of the low-frequency Alfven wa ves. Although the 
observed wind speed around r ~ 10i?s ijGrall et alJl996j) 
appears to exceed our result, it might reflect the wave 
phenomena in the solar wind rather than the outflow 
ilHarmon fc Coles ll2005l) . 

Figure 3 presents the dissipation of the waves. The top 
panel plots the following quantities, 



S c = p5v 



2 (v r + v ph ) 2 r 2 f(r) 



Vph 



(7) 



which is consiste nt with the observed high-speed stream 
around the earth l)Aschwanden. Poland, fc Rabin 1 2001). 



for outgoing Alfven , incoming Alfven , and outgoing slow 
MHD (sound) waves, where Sv and v p h are amplitude and 
phase speed of each wave mode, and S c 's are normalized 
at r c — 1.02. Bs. S r . is an adi abatic constant derived from 
wave action l| Jacques 1119771) . and corresponds to the en- 
ergy flux in static media. S c is conserved in the expand- 
ing atmosphere if the wave does not dissipate, while it 
is not the case for the energy flux. For the incoming 
Alfven wave in Figure 3, we plot the opposite sign of 
S c so that it becomes positive in the sub-Alfvenic re- 
gion. The outgoing and incoming Alfven waves are de- 
composed by correlation between v± and B± . Extraction 
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of the slow wave is also from fluctuating components of v r 
and p. The bottom panel plots (6va,+)/va, ($va,+) /cq, 
and (Svs,+) /cs, where (5va,+) and (5«s,+) are ampli- 
tudes of the outgoing Alfven and slow modes. All the 
quantities in Figure 3 are time-averaged over 30 minutes 
to smooth out the variation due to the phase. 

The top panel shows that the outgoing Alfven waves 
dissipate quite effectively; S c becomes only ~ 10~ 3 of 
the initial value at the outer boundary. First, a sizable 
amount is reflected back downward below the coronal 
base (r — i?s < 0.01i?s), which is clearly illustrated as 
the incoming Alfven wave following the outgoing compo- 
nent with slightly smaller level. This is because the wave 
shape is considerably deformed owing to the steep den- 
sity gradient; a typical variation scale (< 10 5 km) of the 
Alfven speed becomes comparable to or even shorter than 
the wavelength (= 10 4 — 10 6 km). Although the energy 
flux, ~ 5 x 10 5 erg cm~ 2 s _1 , of the outgoing Alfven waves 
(S c in the static region is equivalent with the energy flux) 
which penetrates into the corona is only ~ 15% of the 
input value, it satisfies the requirement for the ener gy 
budget in the coronal holes ijWithbroe fc Noves Ill977|) . 

Second, slow MHD waves are generated in the corona 
(see Sakurai 2002 for observation) as shown in the figure. 
The amplitude of the Alfven waves is amplified through 
the upward propagation. The bottom panel shows that 
(Sva,+)/cs ^ 1 in r > 2i?g, which means that the 
wave pressure (Bj_/8n) exceeds the gas pressure (eq.pj) 
in spite of the weak nonlinearity, (Sv\.+) /v\ ~ 0.1. 
This excites longitudinal slow waves (Kudoh & Shibata 
1999), which eventually become nonlinear (fas -i_) /cs — 1 
at r ~ 3i?s to lead to the shock dissipation l)Suzukil 
2002). The coronal heating and wind acceleration are 
thus achieved by transferring the energy and momentum 
flux of the outgoing Alfven waves. Fast MHD shocks 
by the direct steepening of the linearly polarized compo- 
nents of the Alfven waves also contribute to the heating 
(Hol rwegirmsl ISuzukil 12004^ . though they are less ef- 
ficient because of the weak nonlinearity. The incoming 



Alfven waves are generated in the corona by the reflection 
of the outgoing ones by the density fluctuations due to 
the slow waves. The reflected waves further play a role in 
the wave dissipation by nonlinear wave- wave interaction. 

A key mechanism in the heating and acceleration of the 
plasma is the generation of the slow MHD waves. Thus, 
one of the predictions from our simulation is the existence 
of the longitudinal fluctuations in the solar wind plasma. 
This is directly testable by in situ measurements of fu- 
ture missions, Solar Orbiter and Solar Probe, which will 
approach to ^45 and 4 i?s, respectively, corresponding 
to the inside of our computation domain. 

We have shown by the self-consistent simulation 
that the dissipation of the low-frequency Alfven waves 
through the generation of the compressive waves and 
shocks is one of the solutions for the heating and ac- 
celeration of the plasma in the coronal holes. However, 
the validity of the 1-D MHD approximation we have 
adopted needs to be examined. Generally, the shock 
dissipation tends to be overestimated in the 1-D simu- 
lation because the shocks cannot be diluted by the geo- 
metrical expansion. On the other hand, there are other 
dissipation me chanisms due to the multidimensionality 
ijOfman Il2004[) . such as turbulent casc ade into the trans- 
verse direction l | O ughton et al. 11200 If) and phase mixing 
(Hevvacrts & Pr SM983). Therefore, the waves might 
also be dissipated by paths different from the shocks in 
the real situations. The self-consistent simulations in- 
cluding these various processes remain to be done to give 
the final conclusion. 
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